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Predictions of the extent of self-enrichment in oxygen of 
giant metal-poor H II regions. 
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ABSTRACT 

In general, Hn regions do not show clear signs of self-enrichment in products from 
massive stars (M > 8M©). In order to explore why, I modeled the contamination with 
Wolf-Rayet (WR) star ejecta of metal-poor [Z = 0.001) Hn regions, ionised either 
by a 10 6 -M Q cluster of coeval stars (cluster 1), or a cluster resulting from continu- 
ous star formation at a rate of lMayr" 1 (cluster 2). The clusters have Z = 0.001 
and a Salpeter initial mass function (IMF) from 0.1 to 120 M Q . Independent one di- 
mensional constant density simulations of the emission-line spectra of unenriched H II 
regions were computed at the discrete ages 1, 2, 3, 4, and 5Myr, with the photoion- 
isation code CLOUDY (jFerland et al.| [l998). using as input, radiative and mechan- 
ical stellar feedbacks predicted by the evolutiona ry synthesis code STARBURST99 
(jLeitherer et al l Il999t IVazquez fc Leithererl I2005D . Each H II r egion was placed at 
the ou ter radius of the adiabatically expanding superbubble of iMac Low fc McCravl 
). For models with thermal and ionisation balance time-scales of less than 1 Myr, 



and with oxygen emission-line ratios in agreement with observations, the volume of 
the superbubble and the H II region were uniformly and instantaneously polluted with 
stellar ejecta predicted by STARBURST99. I obtained a maximum oxygen abundance 
enhancement of 0.025 dex, with cluster 1, at 4 Myr. It would be unobservable. 

Key words: Hn regions - ISM: abundances - ISM: bubbles - ISM: evolution - stars: 
early- type. 



1 INTRODUCTION 

Single massive stars , which have masses M > 
8M l|Matteuccil l200ll . § 2.1.3), play a seminal role in 
the chemical enrichment of galaxies and the intergalactic 
medium. Indeed, they produce the bulk of the cv-elements 
(O, Ne, Mg, Si, S, Ar, Ca, an d Ti) and they contr ibute to 
some Fe and Fe-peak elements (jFrancois et alj|2004 , § 4) . In 
addition, at solar metalli cit^, most of them synt hesize all 
elements from B to Rb (|Limongi fe ChiefBl l2007h . Finally, 
at very low metallicity, i.e., at 12 + log(Q/ H) < 7, they 
prod uce signific ant amounts of C and N (see IPettini et al.l 
|200S| for C; see IChiappini et all 120061 for N), although at 
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1 Metallicity, Z, is the mass fraction of metals, i.e., of ele- 
ments heavier than He. Throughout this work, and for con- 
sistency with the STARB URST99 models, the ref erence solar 
metallicity is Z e ~ 0.020 llAnders fe Grev cssc 1989). bmce oxy- 
gen is the most abundant metal in the universe, the number 
density ratio O/H is often used as a gauge of metallicity. For 
Z e ~ 0.020, 12 + log[O/H] ~ 8.90. Note that a more r ecent 
value is 12 + log[O/H] = 8.66 ± 0.05 jAsplund et alj|2004) . 



the rest of metallicities, their contribution to the latter 
two elements, wh ich are essential to life, is still under 
inves tigation ( see ICarigi et al" | _|2005l versus Kub ono et al.l 
I2006J for C; see lHenrv et al.ll2006] for N). 

Unfortunately, t he chemical yields of such stars are 
rather uncertain (e.g-. lMevnet et a"flll994T ). This issue should 
be resolved as we learn more about stellar mass loss from 
theory and observations (e.g., as in iFullerton et ail 120061 ) ■ 
and as we improve our understanding of ho w specific burnin g 
reaction rates affect nucleosynthesis (e.g.. pTur et al.|[2007l ). 
In addition, observations of environments which have be- 
come polluted with massive star ejecta, on space-scales of 
a few hundred parsecs, and over relatively small time-scales 
(t < 10 7 yr), i.e., of self-enriched environments, should con- 
tribute to calibrate massive star yields. 

Targets associated with massive stars which could be 
observed in order to accomplish this task are: the circum- 
stellar nebulae produced by the winds of red supergiants 
(RSGs) or luminous blue variables (LBVs), the circumstel- 
lar bubbles formed by the winds of WR stars, supernova 
remnants (SNRs), the superbubbles formed by the winds 
and supernovae (SNe) of massive stars in associations, and 



2 Aida Wofford 



Hn regions, in particular, low metallicity Hn regions, since 
self-enrichment is expected to be more noticeable when the 
content in metals of the 'to be contaminated' gas is initially 
low. 

Although our focus is on H II regions, which are among 
the most reliable targets for obtaining Galactic and extra- 
galactic chemical abundances, the rest of the above objects 
may affect the spectra of giant and extragalactic H II regions. 
For this reason, let us briefly review enrichment in all of 
them. 

After the main sequence (MS), massive stars be- 
come RSGs, or if they have very larg e birth masses 
(M ~ 60 M ), then they become LBVs (|Nota fc Lamersl 
1997). During these stages, stars turn on slow but dense 
winds. These winds slowly expand away, producing cir- 
cum stellar nebulae enriched w ith processed stellar mate- 
rial (jGarcfa-Segura et ai1ll996al lbT; ISmith et al.lll998l ; iGehr j 
|2004| ; iHirschi et al.1 20071 ). Later o n, if the star's birth mas s 
is M > 30 M at solar metallicity (|Limongi fc Chiem|[2007l ) . 
or higher at a lower metallicity, then the star will become a 
WR star. During this stage, the star turns on an enriched 
wind again. This fast wind sweeps up the previous circum- 
stellar nebula, whi ch collects in a circumstella r s hell known 
as a wind bubble JGarcfa-Segura et al.l fl996al ibl; IChu et all 
l2004lArthurll2007r i. 

Moving to mass loss at the end of the stellar lifetime, 
X-ray observations can show e nrichment in individual young 
SNRs. iRakowski et al.l (|2006l ) used radio and X-ray obser- 
vations of the ejecta-rich Galactic SNR G337.2-0.7 in or- 
der to determine that the origin of the remnant is a SNIa 
and not a core collapse SN. In the optical, although abun- 
dance studies of SNR s and H11 regions in the s piral galax- 
ies M 33 and M31 by iDopita et all (|l980h and iBlair et all 
(1982), respectively, revealed no signs of enrichment in N at 
a given radi us, recent optical obse rvations of SNRs show N- 
enrichment dRudie fc Fesenll2007l ), as well as O-enrichmcnt 
l|Ghavamian et aT 20051 ) from core collapse SNe. 

Finally, when multiple stars of spectral type O and 
B are gravitationally bound in a cluster or association, 
then their individual circumstellar shells expand inside the 
superbubble blown collect i vely by their winds and SNe 
llMac Low fc McCravl Il988l ; IChu et al.1 12004 ICho fc Kangj 
2008). The interiors of such superbubbles are filled with hot 
(10 6 - 10 7 K) gas whose origi n is shocked fast w ind which 
then mixes in with SN ejecta. ISilich et all l|200ll ) predicted 
that superbubble enrichment should be observable within 
the first 10 r yr of the superbubble's lifetime. However, no 
convincing cases of self-enrichment have been reported from 
X-ray observations of the hot gas inside superbubbles in our 
Galaxy or in the Magellanic Clouds. One plausible expla- 
nation is that some of the interstellar gas swept up by the 
superbubble, which has a temperature of T < 10 4 K, i.e., 
which is cooler than the million-degree gas inside the su- 
perbubble, evaporates into the hot superbubble's interior, 
diluting its enriched content. 

We now return to Hn regions. These ionised nebulae 
can be found around circumstellar bubble s and superbub- 
bles l|Townslev et al.ll2003l ; IChu et al.ll2004l ). First we review 
some of the evidence showing that these objects do not be- 
come significantly polluted with massive star ejecta during 
their lifetime (t ~ 10 7 yr). Measurements in the Galactic 
disk, as well as in other spirals, indicate that fluctuations 



in the metal abundances of H 11 regions located at the same 
galactic radius are less t han 20 — 30 per cent on average 
(|Scalo fc Elmegreen|[2003 ). For instance, based on a study 
of 41 Hn regions in MIOl. lKennicutt fc Garnettl (|l996l ) sug- 
gest that the dispersion in the log(0/H) abundance is con- 
siderably less than ±0.1 — ±0.2 dex. In addition, localized 
enrichment does not occur in the Orion nebul a, which is only 
a few hundred parsecs away from us. Indeed. [Esteban et al.l 
(2004) measured heavy element abundances in this Hn re- 
gion and found that they are just slightly abo ve solar. Re- 
gardin g low metallicity Hn regions, a study by IKobulnickvl 
(1999) of nearly 100 low-mass, metal-poor starburst galax- 
ies, aimed at detecting He-, N- and O-pollution by WR stars, 
in Hn regions, revea led no significant variations over spa- 
tial scales I > 10 pc. IKobulnickvl (1999) offered the follow- 
ing three explanations for the lack of success in the study: 
(1) an IMF abnormally truncated at the high-mass end, 
or steeper, which would reduce the amount of contaminat- 
ing ejecta from massive stars, (2) stellar black holes, which 
would swallow the polluting ejecta, (3) ejecta predominantly 
located in the cold phase (T = 10 2 K), or alternatively, in 
the hot phase (T = 10 6 K) of the interstellar medium (ISM). 

IKobulnickvl (|l999l ) rejected the idea of an 'abnormal' 
IMF, on the basis of the undeniable presence of massive 
WR stars. Stellar black holes on the other hand are always 
a possibility. It would be interesting if they can be detected. 
Finally, the hypothesis that the ejecta is predominantly lo- 
cated in the hot phas e of the ISM, is the basis of the model of 
iTenorio- Tagle|[| 19961 ). Indeed, in the latter model, the ejecta 
from OB stars in associations remains trapped within the 
supperbubble carved in the ISM by stellar winds and type II 
supernova (SNII) explosions, while the superbubble expands 
through the ISM. This is because two gases can only mix 
through diffusion efficiently if they have similar high tem- 
peratures and similar densities, but the conditions in the su- 



perbubble interior (T = 10 b - 10 7 K, n = 10" 
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and in the immediate ISM, i.e., in the ISM swept- up by the 
superbubble (T < 10 4 K, n = 10° - 10 2 cm" 3 ), are different. 
The ejecta is only released by the superbubble when the su- 
perbubble reaches the halo of its host galaxy and breaks up 
due to the negative density gradient. The ejecta may then 
blow out of the galaxy or it may rain back in the form of col d 
'droplets' of 10 13 - 10 15 cm in size (|Stasiriska et al.l 120071 ). 
where molecules and grains have formed. These droplets 
spread over kiloparsec scales and mix at the atomic level 
with the ISM when the radiation from a new generation of 
massive stars dissociates the droplet molecules and grains, 
10 8 yr after the matter was ejected from the stars. Note that 
interstellar metal enhancements, such as the droplets, would 
be hard to observe, given their size, in particular if signifi- 
cant dilution of the stellar ejecta with interstellar material 
occurs within the superbubble. 

However, there is some evidence that Hn regions be- 
come polluted with massiv e star ejecta, during the i r life - 
times, at most metallicities. iKobulnickv fc Skillmanl (|l996l ) 
found a fluctuation of 0.095 ± 0.019 dex in log(0/H) 
over a scale of 200 pc that is consistent with pollution 
from SNell, betwe e n two starburst regions of NGC 4214. 
IStasihska fc Izotovl |2003l) detected possible enrichments in 
N and maybe in O in a study of 400 H II galaxies spanning 
the metallicity range Z/Zq = 0.02 — 0.2. The oxygen abun- 
dances recently determined within an area of 400 x 300 pc 2 , 
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in the H II galaxy II Zw 70, sp an the range 12 + log(0/H) = 
7.65 - 8.05 |Kehrig et ai1l2008l) . and the contribution of pol- 
lution with massi ve star ejecta to this s catter is still to be 
determined. Also, iTsamis fc PequignotJ |2005h have argued 
that 30 Doradus in the LMC, which is the closest giant Hn 
region to us, is affected by pollution from SN ejecta. Further- 
more, a recent model of a wind and radiation driven bubble 
around a solar metallicity star of 85 Mq , shows that at the 
end of the stellar lifetime, the amount of C reaches 22.3 per 
cent above solar in the associated H II region, although the 
amounts of O and N, in the same location, are 'insignifi- 
cantly higher than solar' (|Kroger et al ] |2006h . Finally, the 
most well established case of H II region self-enrichment is 
that of the starburst galaxy NGC 5253, where N was re- 
ported to be larger by a factor of 3 at two locations with 
respect to a third one, within a length-scale of 50 — 60 pc 
|Kobulnickv et all 1 19971 ). This enrichment, which seems to 
be accompanied by an enhancement in He, in agreement 
with massive star evolution , was recently re-established by 
lLopez-Sanchez et al.l l|2007t ) , and is probably associated with 
the circumstellar environment of underlying WR stars. Note 
that galaxy Markarian 996 cou ld be a case similar to that of 
NGC 5253 (|Thuan et al.lll996T ). In summary, the extent to 
which H II regions become self-enriched during their lifetime 
is still under investigation. 

As can be seen from the above summary, predictions 
of the physical and spectral properties of self-enriched H II 
regions, which are useful for learning more about the fate 
of massive star products, can only be as good as our un- 
derstanding of how such stars evolve and affect their cir- 
cumstellar environment. In order to move forward in our 
understanding of the chemical enrichment of the ISM, I pre- 
dicted oxygen abundance enhancements, due to pollution 
with massive star ejecta, of metal-poor (Z = 0.001) Hll 
regions, with ages t < 5 Myr. The reason for studying oxy- 
gen is that this element is the dominant metal in massive 
star ejecta and a dominant coolant in Hll regions, there- 
fore, it produces strong observable spectral lines. The reason 
for choosing Z = 0.001 is that self-enrichment is expected 
to be more noticeable when the initial content in metals 
of the ISM is low. The reason for selecting ages less than 
5 Myr is that I concentrated on pollution with WR star 
ejecta. I considered clouds of different gas densitites and 
ionised by star clusters resulting from different star forma- 
tion laws (SFLs). I analyzed the following steady-state prop- 
erties of these clouds: (a) time-scale to reach thermal and 
ionisation balance; (b) values of the [O i] A6300/[O n] A3727 
and [On] A3727/[0 in] AA4959, 5007 emission-line strength 
ratios, R12 and R23, respectively; (c) thickness of the H II re- 
gion; and (d) contribution Alog(0/H) of massive star ejecta 
to the total oxygen abundance of the Hll region. I stud- 
ied how the evolution of the spectral energy distribution of 
the stars, the increase in the separation between the Hll 
region and the stars, and the dust content and hydrogen 
density of the Hll region, affect properties (a) to (d). As a 
test of the validity of each model, the ratios R 12 and R23 
were c ompared with observed values collected bv lNava et alj 
(2006). The paper is organized as follows. In § [2] I outline 
the capabilities of the two computer codes used in this study 
and provide the input parameters and model assumptions. 
In § El I discuss the output from the computer codes and 
explain my results for (a) to (d). In § [4[ I summarize and 



conclude. In order to make th e paper easier to read , detail s 
on the superbubble model of iMac Low fc McCravl l|l988h . 
and part of the discussion on the model assumptions are 
provided in Appendices A and B. 



2 MODELS 

Here I outline the capabilities of the two computer codes 
used in my work and I provide the input parameters and 
model assumptions. 

2.1 STARBURST99 

STARBURST99 llLeitherer et all Il999l ; 

IVazauez fc Leithererl 120051 ). is a widely used evolu- 
tionary synthesis code for modeling various properties 
of stellar populations. Initially, the code was optimized 
for young stellar populations (t < 10 8 yr), which coined 
its name. However, the current running version (5.1) 
allows to model older stell ar populations, as described in 
IVazauez fc Leithererl (|2005l ). Different choices are available 
for the following input parameters: star formation law and 
related quantities, metallicity and stellar evolution tracks, 
and stellar atmospheres and spectra. I used version 5.1 
of the code for predicting the radiative, mechanical, and 
chemical feedbacks of ionising stellar populations, between 
0.1 and 5 Myr. Here is my choice of input parameters. 

2.1.1 Star formation law 

I tested two s tar fo rmation laws, which according to 
iLeitherer et al.l (1 19991 ), reproduce well typical spectropho- 
tometric properties of star-forming regions in galaxies, i.e., 
a single burst of star formation with a total birth mass of 
10 6 Mq , where the stars are born simultaneously and instan- 
taneously (hereafter, instantaneous SFL), and continuous 
star formation at the constant rate of 1 Mq yr _1 (hereafter, 
continuous SFL). For instance, a 10 6 Mo-cluster is thought 
to be responsible for an ultra-compact Hll r egion in star- 
burst galaxy NGC 5253 (iTurner fc Beckll2004l ). Henceforth, 
the stellar populations resulting from the instantaneous and 
continuous SFLs will be referred to as clusters 1 and 2, re- 
spectively. 

2.1.2 Initial mass function 

Stars span the mass range from 0.1 to 120 Mq and ar e dis- 
tributed according to a Salpeter IMF (Salpeter 1955), i.e., 
such that the initial number of stars per unit volume (in pc 3 ) 
and with masses between m and m + dm (in Mq) is given 
by dN = k m -2,35 dm, where the constant k, normalizes the 
mass distribution to the total stellar mass. Given this IMF, 
the initial number of O V stars in cluster 1 is ~2xl0 3 . 

2.1.3 Final stages of stellar evolution 

The minimum mass of a SNII progenitor was set to 8 Mq , 
and the cut-off mass for the formation of black hole to 
120 Mq, i.e., no black holes are formed. The choice of no 
black holes maximizes the total amount of contaminating 
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material for the H II region, since the ejecta from supernovae 
are not swallowed by the remnants. 



2.1.4 Metallicity and tracks 

I used the set of evolutionary tracks best suited for young 
stellar populations, i.e., the set of non-rotating Geneva 
tracks. Since self-enrichment is expected to be more no- 
ticeable at low metallicity, I selected the lowest metallicity 
available for these tracks, i.e., Z = 0.001. The adopted set of 
tracks uses the standard mass loss models of ISchaller et al.l 
<| 19921 ) for stars with mass es M < 12Mq, and the high 
mass loss models of iMevnet et al.l (Il994h for s tars with 



Mevnet et all (|l994h cover all 



M > 12 M Q . The tracks of 
phases up to the end of central C burning and reproduce 
well most of the observed properties of WR stars, in partic- 
ular, the WR subtype distributions and the low luminosities 
reached by certain WR stars, alt hough they do not rep ro- 



r 



duce well the mass loss rate itself l)Leitherer et al.ll 199 



2.1.5 Chemical yields 

The wind chemical yields, i.e., the mass in each element 
ejected by winds in M yr _1 , come straight out of the evo- 
lutionary tracks. Note that since only H, He, C, N, and 
O have significant changes in the pre-supernova phase, the 
wind yields for elements like Si or Fe are simply determined 
by the initial abundan ces. On the other hand, the yields 
from SNell are those oflWooslev fc Weaver! (|l995T ). with the 
corrections of lLindner et al-H 1999T ) . 



2.1.6 Stellar atmospheres and spectra 

In STARBURST99, the ensemble of stars is treated as a 
single source with a luminosity equal to the sum of the lu- 
minosities of the individual stars, and with a spectral shape 
equal to the superposition of the spe ctra of the sta r s. The 
individual stellar spectra come from iLeieune et alj (Il997 ) 
for st ars with plane-p arallel atmospheres, Pauldrach et al.l 
i|200ll ) for O stars, and lHillier fc Milled (|l998l ) for other stars 
with strong winds. The cluster spectra span the wavelength 
range 91 — 1.6 x 10 6 A. They are composed of 1220 points, 
which are spaced such that AA = 2 — 10 A for A < 912 A 
and AA = 10 - 2 x 10 5 A for A > 912 A. 



2.2 CLOUDY 



CLOUDY (jFerland et all Il998h . is a widely used 
plasma/molecular photoionisation code that predicts 
the thermal, ionisation, and chemical structures of clouds, 
ranging from intergalactic to the high-density local- 
thermodynamic-equilibrium limit, and their spectra. The 
code requires the following input parameters: 1) shape 
and luminosity of the external radiation field striking the 
cloud (hereafter, the incident continuum), 2) chemical 
composition and grain content of the gas, and 3) geometry 
of the gas, including the radial extent of the gas and the 
dependence of its density on the radius. I computed one 
dimensional constant density steady-state models of the 
emission-line spectra of Hn regions, with version 07.02 of 
CLOUDY, using for 1) the spectral energy distributions 



predicted by STARBURST99 at times t = 1, 2, 3, 4, and 
5Myr since the beginning of star formation, and for 3) 
the mechanical feedback of the clusters, as explained in 
§ 12.2.21 Note that the H 11 region models are independent 
of each other, i.e., the steady-state properties of a given 
cloud do not affect the properties of subsequent models. 
The computation of time-dependent H 11 region simulations 
may be possible with the 2009 release of CLOUDY (G. 
Ferland, private communication). Further model details are 
provided below. 

2.2.1 Gas and dust chemical compositions 

The Hn region gas is composed of the 15 most abundant 
elements by number in the sun. All of these elements, ex- 
cept for Na, Al and Ca, have lines in the wavelength range 
3, 100 —10, 400 A that have been observed in Galactic H 11 re- 
gions l|Garcfa-Roias et al.ll2004h . On the other hand, Na, Al 
and Ca are refractory ele ments expected to be present in the 
form of dust in the ISM (|Sembachll 19951 : ISembach fc Savage! 
ll996l : lFerlandll2006h . The exclusion of additional elements 
should have little impact on the results, given the low metal- 
licity of the clouds (Z ~ 0.001). The abundances of ele- 
ments in the gas and dust, used as input to the CLOUDY 
models, are given in columns 7 and 8 of Table [TJ respec- 
tively. They were computed as follows. The gas abundances 
were obtained by applying the scaling factors in column 5 to 
the solar values in column 2. Similarly, the dust abundances 
were obtained by applying the scaling factors in column 6 to 
the values in column 3, which are CLOUDY's default abun- 
dances for the Galactic ISM. Column 4 of Table [T] shows 
that C and O are present in CLOUDY's dust, and that dust 
has a high content of C. However, because dust abundances 
are highly uncertain, I did not deplete C and O from the 
gas in order to account for their presence in the dust, as I 
did for less abundant elements. In § [3] I show that the ef- 
fect of dust on the emitted spectrum of the H 11 regions in 
negligible. This is due to the low metallicity of the models. 

For consistency between CLOUDY and STAR- 
BURST99, the helium mass fraction Y was computed from 



Y = Y P + 



Yq-Yp 
Zq — Zp 



(1) 



where the primordial values are Yp = 0.24 and Zp — 0, 
while the solar valu es are Yq — 0.299 and Z — 0.019 
(|Schaller et all 1 1992ft . Equation (TJ yields Y = 0.243. The 
ratio of the helium to hydrogen number abundances, y = 
n(He)/n(H), was derived from Eqn. (flj and 



Y = 



M(He) 



Ay 



M (H) + M(He) + M(Z) " 1 - Ay ' ^ 

where M(H), M(He), and M(Z) are the masses in H, He 
and metals, and the right hand side of the Eqn. ([2} comes 
from neglecting M(Z) and using M(H) = m(H) n(H) V and 
M(He) = 4m(H) n(He) V, where m(H) is the atomic mass 
of H, and V is some volume. The resulting He abundance is 
n(He) » 0.1n(H). 



2.2.2 Geometry 

The ionising cluster is located at the center. The cloud has 
a covering factor of ^An (i.e., the central source is fully 



Predictions of the extent of self- enrichment in oxygen of giant metal-poor Hn regions. 5 

Table 1. Gas and dust abundances in units of 12 + log[n/n(H)], and scaling factors applied to n/n(H). 



element 




reference value 


scaling 


factor 


value of model cloud 


sun 2 


dust 


dust / sun^ 


gas 5 


dust 6 
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■ 






12.00 




He 7 


10 99 










10.90 




U 


O.OO 


8.44 


0.76 


U.UO 




7.26 


7.14 


N 


8.05 






0.05 




6.75 







8.93 


8.12 


0.15 


0.05 




7.63 


6.82 


Ne 


8.09 






0.05 




6.79 




Na 


6.32 






0.05x0.1 




4.02 




Mg 


7.58 


7.52 


0.87 


0.05x0.1 


0.05 


5.28 


6.21 


Al 


6.47 






0.05x0.1 




4.17 




Si 


7.55 


7.52 


0.93 


0.05x0.1 


0.05 


5.25 


6.21 


S 


7.24 






0.05 




5.94 




Ar 


6.56 






0.05 




5.26 




Ca 


6.35 






0.05x0.1 




4.05 




Fe 


7.59 


7.52 


0.85 


0.05x0.1 


0.05 


5.29 


6.21 


Ni 


6.25 






0.05x0.1 




3.95 






Figure 1. Model geometry. Hashed shell: bounding shell of su- 
pcrbubble. r^: inner radius of superbubblc's shell, ro: outer radius 
of superbubble's shell. Region between ro and r^'- model cloud. 
Region between rrj and ri: Hn region. Dashed-circle: Hn region's 
ionisation front. 

covered), a filling factor of one (i.e., there are no clumps in 
the gas), a constant value of n(H) (cm -3 ), and it starts at 
distance ro from the center. The value of rp increases w ith 
time as given by the model o f lMac Low fe McCravl[|l988h for 
the adiabatic expansion of a superbubble [see Appendix A 
for a summary]. The outer radius of the superbubble, ro, is 
given by 

r = (267 pc) (3) 

where L38 (10 38 erg s -1 ) is the average mechanical luminos- 
ity from the stellar winds and supernovae (I used the average 
from 0.1 to 5Myr), tj (10 7 yr) is the time since the onset of 
stellar mass loss (it starts on the main sequence for massive 
stars), and no (cm -3 ) is the density of the ambient undis- 
turbed ISM at the location of the superbubble (I computed 
models with no = 1 and no = 10). Note that by letting 
the Hn region start at ro, I am assuming that its ionisation 
front is located outside of the superbubble. In nature, the 
ionisation front may be located within the superbubble. In 



§ [21 I compare the value of ro with the Stromgren radius of 
each Hn region. 

The radial extent of the cloud is set by the point where 
the electron temperature of the gas has dropped to 100 K. 
For most models, the resulting geometry of the emission 
is plane-parallel. Figure [1] is a sketch of the gas geometry. 
As shown in the figure, the thickness of the superbubble's 
bounding shell is ro — rj, the thickness of the model cloud is 
Ti — ro, and the thickness of the Hn region is ri — ro. 



2. 2. 3 Volume of radius ro 

In the superbubble model of iMac Low fe McCravl i|l98St ). 
the volume of radius ri contains gas with a temperature T 
given by 

T = T c (l-r/r ) 2/5 , (4) 
where 

T c = (3.5 x 10 6 K)L 8 g 35 n /35 t~ 6/35 (5) 

is the central temperature; and an atomic density n given 
by 

n = n c (l-r/r )" 2/5 , (6) 
where 

I A „ in-3 -3\ r 6/35 19/35 .-22/35 /~\ 

n c = (4 x 10 cm )L 3 ' S n t 7 (7) 

is the central density. I assumed that the value of n at rj is 
constant for n < r < r2, and I used Eqns. Q to (J7J) to derive 
the hydrogen density n(H) of the model cloud, as explained 
next. Note that for the purpose of radiative transfer, the 
volume of radius ro, Vo, was treated as a cavity empty of 
gas. This assumption is discussed in Appendix B. 



2.2.4 Hydrogen density 

For simplicity, I assumed that n(H) is constant throughout 
the model cloud. However, I assumed that n(H) varies with 
time as given by 
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»(H) = 



1.4 x L 
1.1 



2/5 3/5 



'0 



,4/5 
t 7 



(«) 



I obtained Eqn. ([SJ as follows. First I solved Eqn. Q for 
the radius, which yields, r = r [l - (T/T c f /2 \ . Then I 
substituted the latter expression in Eqn. (JSJ), which yields 
n = n c T c /T. Next, I used n ft! n(H) + n(He) fa 1.1 x n(H). 
Finally, I assumed that T = 10 4 K, which is the characteris- 
tic equilibrium temperature of an Hn region. In summary, 
the atomic density of the model cloud was assumed to be 
equal to the atomic density given by Eqn. (|6}, at the point 
where T = 10 4 K. 



2.3 Self-enrichment 

Let O/H = n(0)/n(H). In addition, let i and / stand for 
'initial' and 'final', where initial is 'prior to' and final is 'af- 
ter' contamination of the cloud with stellar ejecta. Finally, 
let An be the increase in the number density of an element 
in the Hn region, due to contamination with massive star 
ejecta. The contribution Alog(0/H) of the ejecta to the oxy- 
gen abundance of the H II region, is given by 



where 

O 
H 



m(Q) + An(O) 
rii(H) + An(H) ' 



(9) 



(10) 



The value of rn(H) is given by Eqn. JS), that of rii(O) by 
column 7 of Tabled] and that of An(X) by 



An(X) = 



M(X) 



(Vb + Vi)m(X)' 



(11) 



where M(X) is the cumulative mass in element X ejected 
by the stars, up to the time under consideration, Vo is the 
volume of radius ro , Vi is the volume of the shell of thickness 
n — ro, and m(X) is the atomic mass of element X. 

At this point, one may wonder what is the physical 
mechanism assumed for the mixing of the stellar ejecta 
with the Hn region. Turbulent mixing is expected to oc- 
cur in shock fronts and contact discontinuities, where gas 
is stretched and folded until the gradient length appro aches 
the collision mean free path |Scalo fc Elmegreenl feOCM') . I as- 
sumed that stellar winds and supernovae disperse the ejecta 
throughout the volume of radius fj, and that turbulence in 
the wind and SN shocks mixes the ejecta with the gas pre- 
existing in the volu me, quickly and effic iently. On the other 
hand, according to iTenorio-Tagiel { 1996), mixing across the 
superbubble's shell is inefficient. This is because two gases 
can only mix through diffusion efficiently if they have similar 
high temperatures and similar densities, but the conditions 
of the gas for r < n (T ~ 10 6 - 10 7 K, n ~ 10" 3 - 1(T 2 
cm -3 ) and in theHu region (T ps 10 4 K, n = 10-10 2 cm -3 ) 
are different. Without providing a mechanism for the mixing 
across the superbubble's shell, I test if uniform and instanta- 
neous mixing of the stellar ejecta, with the gas in the volume 
Vo + Vi, results in an observable enhancement in the O/H 
abundance of the Hn region. 




Figure 2. Number of O stars [iV(O), panel a] and WR stars 
[AT(WR), panel b] versus time (t). Curves: solid, cluster 1; dashed, 
cluster 2. Time-step: 0.1 Myr. 



2.3.1 Supplementary models 

For each cluster and value of no, three supplementary models 
were computed at 4 Myr: one without dust, one where the 
hydrogen density was arbitrarily set to n(H) = 100 cm" 3 , 
and one where the distance to the cloud was arbitrarily set 
to ro = 100 pc. The purpose of these models is to test the 
sensitivity of Alog(0/H) to the dust content of the cloud, 
n(H), and ro. 



3 RESULTS 

Here I discuss properties of the radiative, mechanical, and 
chemical feedbacks of clusters 1 (instantaneous SFL) and 
2 (continuous SFL). Hereafter, superindexes 1 and 2 are 
used for quantities related to each cluster, respectively. In 
this section, I also describe the output from CLOUDY, 
compare it with observations, and discuss the contribution 
Alog(0/H) of massive star ejecta to the total oxygen abun- 
dance of the Hn regions. 



3.1 Number of O and WR stars 

Figure [2] shows evolutionary plots of the number of O stars 
[N(0), panel a] and WR stars [JV(WR), panel b], for clus- 
ters 1 (solid curve) and 2 (dashed curve). As time increases, 
iV(O) 1 decreases while JV(O) 2 increases. Also, iV(WR) 1 / 
from 3 to 3.7 Myr, and iV(WR) 2 / for t > 3 Myr. Type O 
and WR stars affect the radiative, mechanical and chemical 
feedbacks of the clusters, as described below. 



3.2 Stellar radiative feedback 

Figure [3] shows plots of the monochromatic luminosity (-La) 
versus wavelength (A), for clusters 1 (panel a) and 2 (panel 
b), at the discrete times 1, 2, 3, 4, and 5 Myr, and in the 
region below the Lyman limit (A < 912 A), which is the 
stellar spectral region of interest, since in Hn regions, most 
of the Hi is the ground state, and therefore, the number of 
hydrogen i onisations by starlight photo ns with A > 912 A is 
negligible (|Osterbrock fc Ferlandl fcooe). Let Li be the ion- 
ising luminosity below the Lyman limit. The luminosity 
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Table 2. Radiative, mechanical, and chemical feedbacks of clusters 1 and 2, at ages 1-5 Myr, predicted by STARBURST99. Entries (all 
logarithmic values except for time): time since the first stars formed, t; number of ionising photons in the Hi, He I, and Hell continua, 
Q(Hl), Q(Hel), and Q(Hcll); bolometric luminosity, Lboli mechanical luminosity from winds and supernovae, Z/ mec h; an d masses in 
hydrogen and oxygen ejected by the stars up to time t, M(H) and M(0). 



t 


Q(Hi) 


Q(Hei) 


Q(Hcn) 


Lbol 


Lmech 


M(H) 


M(O) 


Myr 


s" 1 


s- 1 


s" 1 


erg s 


erg s 


M 


M 


cluster 1 


1 


52.563 


52.034 


48.321 


42.398 


38.467 


1.934 


-1.222 


2 


52.621 


52.093 


48.570 


42.454 


38.980 


2.505 


-0.638 


3 


52.499 


51.843 


47.907 


42.527 


39.392 


3.301 


-0.009 


1 


52.224 


51.463 


47.368 


42.364 


39.983 


3.732 


3.462 


5 


51.959 


51.056 


46.224 


42.203 


39.989 


3.833 


3.462 


cluster 2 


1 


52.579 


52.051 


48.341 


42.415 


38.375 


1.602 


-1.553 


2 


52.888 


52.359 


48.672 


42.723 


38.900 


2.342 


-0.824 


3 


53.061 


52.511 


48.811 


42.923 


39.397 


2.987 


-0.187 


4 


53.142 


52.574 


49.494 


43.050 


40.390 


3.672 


3.176 


5 


53.179 


52.595 


49.495 


43.117 


40.537 


4.041 


3.643 



-r~ 38 



-T* 34 



1 














1 Myr 






~- A 




_ a 


2 Myr 

3 Myr 

4 Myr 

5 Myr 












/ / 
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Figure 3. Monochromatic luminosity (L\) versus wavelength 
(A). Panels: a, cluster 1; b: cluster 2. Curves: thick-solid, 1 Myr; 
dashed, 2 Myr; dotted-dashed, 3 Myr; thick-dashed, 4 Myr; thin- 
solid, 5 Myr. Vertical-dotted lines at 228, 504, and 912 A, ap- 
proximate ground-state ionising wavelengths of Hen, He I and 
Hi, respectively. 




38 - 



1x10 S 2x10 6 3x10 S 4x10 S 5x10 6 

t(yr) 

Figure 4. Mechanical luminosity from stellar winds and SNc 
(Lmech) versus time (t). Panels: a, cluster 1; b, cluster 2. Curves: 
thick-solid, winds and SNe; thin-solid, OB stars; dashed, WR 
stars; dotted-dashed, SNe (lb and II). Horizontal-dotted lines: 
average values log (£ me ch)i from 0.1 to 5 Myr, for clusters 1 and 
2. Time-step: 0.1 Myr. 



is dominated by the most massive O stars during the first 
Myr, by WR stars from 3 to 3.7 Myr, and by less massive 
O stars from 3.7 to 5 Myr. The value of L\ decreases with 
time as the stars contributing to it fade away. On the other 
hand, Lf increases with time for the first Myr, and then be- 
comes independent of time, when a balance is established 
between the birth and death of the stars contributing to it 
(this is why in panel b of Fig. [3] the thick-dashed and thin- 
solid curves, which correspond to 4 and 5 Myr respectively, 
overlap). Finally, L\ < L\ from 1 to 5 Myr. This is due to 
the continuous formation of ionising stars in cluster 2. In 
particular, the Hell continuum is a lot more prominent for 
cluster 2, from 4 to 5 Myr, due to the strong radiation in 
that spectral region, from WR stars, which are absent at 
those times, in cluster 1. 

Table [2] complements Figs. [2] and [3] The entries in the 
table are: time since the formation of the first stars (t), num- 
ber of ionising photons in the Hi, He I, and Hell continua, 
[Q(Hl), Q(Hel), and Q(Hell) respectively], bolometric lu- 



minosity (Lboi), mechanical luminosity from winds and SNe 
(Lmech), and mass in oxygen ejected by the cluster up to the 
time considered [Af(O)]. The evolutions of L moc h and AZ(O) 
are discussed next. 

3.3 Stellar mechanical feedback 

Let L mec h be the mechanical luminosity from stellar winds 
and SNe. Values of L rnech and L mech at 1, 2, 3, 4, and 5 Myr 
are given in columns 6 and 12 of Table[2j respectively. These 
columns show that L mech > L mech at 1 and 2 Myr, whereas 
Lmoch < L moch at 3, 4, and 5 Myr. Figure |4] shows plots 
of L mcc h versus time, for clusters 1 (panel a) and 2 (panel 
b). The contributions from different components have been 
plotted separately. The luminosity L mech is dominated by 
winds from OB stars from 1 to 3 Myr, by winds from WR 
stars from 3 to 3.7 Myr, and by SNe of type lb and II from 4 
to 5 Myr. On the other hand, L mcc h is dominated by winds 
from OB stars from 1 to 3 Myr, by winds from WR stars 
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Figure 5. Oxygen mass loss rate [dM(0)/dt] versus time (t). 
Curves: solid, cluster 1; dashed, cluster 2. Time-step: 0.1 Myr. 



from 3 to 4.8 Myr, and by SNe from 4.8 to 5 Myr. The 
horizontal-dotted lines in Fig. [4] indicate the average val- 
ues of log(L mec h), from 0.1 to 5 Myr, for clusters 1 and 2. 
Recall that the reason for adopting the average (1/mech), 
is that Eqn. ((3} for the radius ro is based on a constant 
input of mechanical energy from the stars. The values of 
^38 = (imcch) /10 38 are Lj 8 « 68 and L| 8 » 97. A com- 
parison of the thick-solid curve and horizontal-dotted line 
in panel a of Fig. [4] shows that L^g is overestimated for 
t < 3 Myr and underestimated for t > 3 Myr. The same is 
true for L§ g , except that the limiting time is t 3.2 Myr. This 
means that ro is underestimated and Alog(0/H) is overes- 
timated at times when WR stars dominate the enrichment 
of the ISM. As shown in 5 13.5.61 that Alog(0/H) is overes- 
timated at the latter times does not change the main con- 
clusion of this study. 



3.4 Stellar chemical feedback 

Figure [5] shows plots of the oxygen mass loss rate dM(0)/dt 
( Mq yr _1 ), for clusters 1 (solid curve) and 2 (dashed curve). 
The value of [dM(0)/dt] 1 slowly increases from 1-3 Myr 
when the winds from OB stars dominate the mass loss. Then 
[dM(0)/dt] a rapidly increases when WR stars appear, and 
it maintains high values while these stars are present (from 3 
to 3.7 Myr). Finally [dM^/dt] 1 falls i nto a dip starting at 
t ~ 4. 2 Myr , which is due to the f act that lWooslev fc Weaver! 
i|l995h and iLindner et all (|l999h did not publish SN yields 
for stellar masses higher than 40 Mq , or equivalently, stel- 
lar ages below 5.6 Myr. The justifications for the dip are 
that SNe with such high masses are very rare and do not 
contribute to the ISM enrichment; and also, that we do 
not know the yields of such events from observations. Note 
that in STARBURST99 v. 5.1, stars more massive than 
40 Mq explode but without any nucleosynthesis, which ex- 
plains why there is not a similar dip in the plot of I^cch • On 
the other hand, [dM(0)/dt] 2 evolves similar to cluster 1 for 
the first Myr, but reaches a steady value around t = 4 Myr, 
which is maintained by the constant presence of WR stars. 



3.5 H II region models 

For each cluster and value of no , Table [3] gives input and 
output properties of the H II region models, as well as some 
supplementary data. The entries in the table are described 
in the caption and are discussed in sections 13.5.11 to 13.5.61 
except for the data corresponding to the 12 supplementary 
simulations at 4 Myr, which either have no dust (models a), 
have n(H) = 100 (models b), or have ro = 100 pc (models 
c). Models a, b, and c are discussed in 



3.5.1 Radius ro 

In order to see if n > ro is a reasonable as- 
sumption, I computed the Stromgren radius, rs = 
{3Q(H)/[4tto4(H,T)]} 1/3 , of each Hn region, where 
Q^(H, T) is the case A radiative recombination coefficient 
to form Hi in any energy level. I used the value at 10 4 K 
(cha racteristic of Hn region s), i.e., a r A = 4.2 x 10~ 13 cm 3 
s~ (lVerner fc Ferland| [i996). For most models, rs > ro, in 
agreement with the assumption that the ionisation front of 
the H II region is located outside of the superbubble. How- 
ever, models 8, 9, 9a, 9b, and 10 have rs < ro. Therefore, the 
latter models are excluded from further analysis. The radius 
ro increases as t increases, no decreases, and L38 increases. 
Continuous star formation yields larger values of ro by a few 
tenths of parsecs, because L§ g > L\ 6 . All values of ro are 
within ~ 100-400 pc, therefore, some of them are compara- 
ble to the thickn ess of the Galactic disk, which is ~200 pc 
|Matteuccill200lh . A superbubble breaks when it reaches the 
halo of its host galaxy, as mentioned in § [T] However, in the 
present work, I assumed that the superbubble and its asso- 
ciated H 11 region remain within the gas-rich component of 
the host galaxy (which is not necessarily a disk galaxy), and 
that stellar ejecta is confined within radius n, from 1-5 Myr. 



3.5.2 Hydrogen density 

The hydrogen density n(H) increases as t decreases, no in- 
creases, and 1/38 increases. All values of n(H) are within 
~10-200 cm -3 , which is in good agreement with observed 
values in giant H 11 regions (e.g., see column 2 of table 5 in 
iNava et alj|2006h . 



3.5.3 Equilibrium, time-scale 

The reason for checking the value of t oq is that in order to 
be valid, simulations must have i eq < 1 Myr. This is be- 
cause subsequent models are 1 Myr apart. For all models, 
the longest process to reach equilibrium is heating-cooling 
balance. Cooling is dominated by collisional cooling due to 
hydrogen lines, which becomes less important as time in- 
creases and n(H) decreases. This is why i cq increases as time 
increases and n(H) decreases. The requirement that i eq be 
less than 1 Myr rules out model 5 as valid. 



3.5.4 Oxygen emission-line ratios 

Figure [6] shows plots of observed versus model values of the 
ratios R12 = 6300/3727 and R23 = 3727/(4959 + 5007). It 
includes the 32 models in Table [3] The filled s ymbols cor- 
respond to the observational data collected by iNava et al.l 
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Table 3. Input and output of the CLOUDY models for the instantaneous and continuous SFLs (cluster 1 and 2, respectively), and 
for the two values of the atomic density no (cm -3 ). The Stromgren radius and value of Alog(0/H) are also provided. Entries: model 
number; age of cloud, t; dust content (0=dust-free, l=dusty); Stromgren radius, rg; outer radius of superbubble, ro; thickness of the 
Hn region, r\ — ro; hydrogen density of the model cloud, n(H); time-scale to reach thermal and ionisation balance, t eq ; emission-line 
strength ratios [O i] A6300/[O n] A3727 and [On] A3727/[0 in] AA4959, 5007; and contribution of the cluster's ejecta to the total oxygen 
abundance of the Hll region, Alog(0/H). 



model 


t 

Myr 


dust 


pc 


ro 
pc 


ri - ro n(H) 
pc cm -3 


tcq 

Myr 


6300 
3727 


3727 
4959+5007 


Alog(0/H) 
dex 












cluster 1, no 


= 1 










1 
2 
3 
4 

4a 
4b 
lc 

5 


1 

2 
3 
4 
4 
4 
1 
5 


1 
1 
1 
1 

1 
1 
1 


891 
932 
849 
687 
687 
687 
687 
561 


156 
236 
301 
358 
358 
358 
100 
110 


10 
15 
13 
8 
8 


62 
5 


43 
25 
18 
14 
14 
100 
14 
12 


0.3 
0.5 
0.7 
0.9 
0.5 
0.1 
0.9 
1.1 


0.02 
0.02 
0.01 
0.01 
0.01 
0.02 
0.01 
0.01 


0.05 
0.06 
0.13 
0.44 
0.44 
7.05 
0.06 
1.91 


0.000 
0.000 
0.000 
0.025 
0.025 

0.227 












cluster 1, no 


= 10 










6 
7 
8 
9 

9a 
9b 
9c 
10 


1 

2 
3 
4 
4 
4 
i 
5 


1 
1 
1 
1 



1 
1 
1 


192 
201 
183 
118 
148 
148 
118 
121 


98 
149 
190 
226 
226 
226 
100 
258 


2 
3 

1 
1 

6 
1 


173 
99 
72 
57 
57 
100 
57 
48 


0.1 
0.1 
0.2 
0.2 
0.1 
0.1 
0.2 
0.3 


0.02 
0.02 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 


0.08 
0.09 
0.21 
0.78 
0.78 
1.66 
0.13 
3.65 


0.000 
0.000 

0.206 












cluster 2, no 


= 1 










11 

12 

13 

14 

14a 

14b 

14c 

15 


1 

2 
3 
4 
4 
4 
i 
5 


1 
1 
1 
1 



1 
1 
1 


902 
1,144 
1,306 
1,390 
1,390 
1,390 
1,390 
1,430 


167 
254 
324 
385 
385 
385 
100 
110 


6 
16 
26 
36 
37 

1 
149 
43 


50 
29 
21 
16 
16 
100 
16 
11 


0.2 
0.4 
0.6 
0.7 
0.4 
0.1 
0.7 
0.8 


0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.03 
0.02 


0.07 
0.05 
0.04 
0.04 
0.04 
0.22 
0.01 
0.04 


0.000 
0.000 
0.000 
0.008 
0.008 
0.002 
0.035 
0.018 












cluster 2, no 


= 10 










16 

17 

18 

19 

19a 

19b 

19c 

20 


1 

2 
3 
4 
4 
4 
4 
5 


1 
1 
1 
1 



1 
1 
1 


194 
246 
281 
299 
299 
299 
299 
308 


106 
160 
204 
243 
243 
243 
100 
277 


1 
3 
1 
6 
6 
3 
26 
7 


199 
114 
83 
66 
66 
100 
66 
55 


0.1 
0.1 
0.2 
0.2 
0.1 
0.1 
0.2 
0.2 


0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 
0.02 


0.11 
0.07 
0.06 
0.06 
0.06 
0.08 
0.01 
0.06 


0.000 
0.000 
0.000 
0.009 
0.009 
0.006 
0.066 
0.021 



(2006), the open symbols to clouds ionised by cluster 1, and 
the crosses to clouds ionised by cluster 2. All models are lo- 
cated on the vertical-dotted lines at 12 + log(0/H) = 7.629, 
which is the input oxygen abundance of the CLOUDY mod- 
els. The value of no is indicated in each panel. The arrows 
shown in panels b and d indicate the direction of time for the 
ten models 1 to 5 and 6 to 10, respectively. For clouds ionised 
by cluster 1, R23 increases as time increases because as the 
most luminous stars in the cluster fade away, the amount of 
+2 decreases, whereas the amount of + increases. 

The observed values in panels a and c present no trend 
and cannot be used to test the validity of the models. Note 
that in most cases, only an upper limit is available for the 
observed value of [O I] A6300. On the other hand, panels b 
and d show that the range of observed values of R23 in- 
creases with the oxygen abundance. This trend was used 



to test the validity of the models. The observational data 
in panels b and d includes no points with R23 greater than 
about 0.45 at 12 + log(0/H) = 7.629. Of the remaining good 
models (exclude models 5, 8, 9, 9a, 9b, and 10), only 4b has 
R23 > 0.45. Therefore, model 4b was also excluded. Note 
that the point with the error bars in Fig. [6] is KISSR 1013 
l|Melbourne et all 12004 ) . whic h is a giant emis sion nebula 
with 12 + log(0/H) = 7.632 (|Nava et alj|2006ft . The hori- 
zontal error bar was computed using standard error propa- 
gation theory and was taken as the characteristic error bar 
of model 4, which is the open circle closest to KISSR 1013 
in Fig. [6] Therefore, model 4 was considered to be a valid 
model. 
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Figure 6. Observed versus model oxygen emission-line strength ratios 6300/3727 and 3727/(4959 + 5007) (including the 32 models in 
Table [3J. Panels: a & b, no = 1; c & d, no = 10. Symbols: filled circles, observed values collected by Nava et al. (2006); open circles, 
clouds ionised by cluster 1; crosses, clouds ionised by cluster 2. Vertical dotted lines: oxygen abundance of model Hn regions. Arrow in 
panel b: direction of time for the five models 1 to 5. Arrow in panel d: direction of time for the five models 6 to 10. 



3.5.5 Thickness of Hll regions 

The radius r\ is given by the point where the fraction of 
ionised to total hydrogen is H + /H = 0.9. The H II regions are 
thin shells with thicknesses (r\ — ro) equal to a few per cent 
of ro. The H II regions ionised by cluster 1 become thinner as 
the ionising stars fade away, while those ionised by cluster 
2 become thicker, as the value of Q(H) increases. 



Of the latter five models, model 4, which corresponds to the 
instantaneous SFL, has the maximum value of Alog(0/H), 
which is 0.025 dex. This value is within the range of obser- 
vational uncertainties i n O/H of metal-p oor H n regions (see 
column 3 of table 6 in lNava et al.|[200q ). In addition, recall 
that the value of ro is underestimated at 4Myr (see § I3.3f> . 
therefore, Alog(0/H) may be even lower. 



3.5.6 Self- enrichment 

The values of An(H) and An(O) were computed from Eqn. 

, using the values in columns 6 and 7 of Table [2j re- 
spectively. For all models, An(H) < 5 x 10~ 3 and n;(H) = 
n(H) > 10. Therefore, Eqn. (|ll|l can be simplified to 



(h)^(h)/ 

which yields 
Alog(0/H) = log 



An(O) 



(h), 



+ 



An(O) 



n(H) 



+ 4.371, 



(12) 



(13) 



Table [2] shows that M(O) is only significant at 4 and 
5 Myr. The fraction of M(O) used to pollute the H II regions, 
is given by Vi/(Va + Vi) = 1 - r'^/rf. It is less than 0.25, for 
valid models, and excluding models c, which have an arbi- 
trary value of ro = 100 pc. Of the models corresponding to 
4 or 5 Myr, models 4, 14, 15, 19, and 20 are good models. 



3.6 Supplementary models 

I proceed to discuss the 12 supplementary models, a, b, 
and c at 4 Myr. The comparisons of models 4 with 4a, 14 
with 14a, and 19 with 19a, show that dust-free clouds (mod- 
els a) are indistinguishable from dusty models. This is due 
to the low metallicity of the clouds. The comparisons of 
models 14 with 14b and 19 with 19b, show that increas- 
ing the value of n(H) by a few xlOcm -3 decreases the 
value of Alog(0/H) by a few xl0~ 3 dex, which is a neg- 
ligible amount. Finally, comparisons of models 4 with 4c, 
14 with 14c, and 19 with 19c, show that reducing ro by 
a few xlOOpc, significantly increases Alog(0/H). In par- 
ticular, the value of Alog(0/H) for model 4 is larger than 
observed oxygen abundance fluctuations in Hll regions of 
higher metalli cities, such as those found in M101 (equal 
to 0.1-0.2 dex iKennicutt fc Garnettl [l99?3) and NGC 4214 
(equal to 0.1 dex iKobulnickv fc Skillmanl fl996l ). Note that 
Eqn. ([3)l for ro is based on a model that neglects the effect 
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of gravity. If gravity was included, this would reduce the 
value of ro, and increase the value of Alog(0/H). 



4 SUMMARY AND CONCLUSION 

I used versio n 5.1 of the evolu t ionar y synthesis code STAR- 
BUR ST99 (|Leitherer et al.1 1 19991 ; IVazquez fc Leithererl 
2005) for predicting the radiative, mechanical, and chemical 
feedbacks of young (t < 5 Myr) metal-poor (Z — 0.001) 
stellar populations, resulting from either a single burst of 
star formation with a total stellar birth mass of 10 6 M 
(cluster 1), or from continuous star formation at a rate 
of IMoyr -1 (cluster 2). I adopted a Salpeter IMF from 
0.1 - 120 M , SNII and black hole cut-off masses of 8 and 
120 Mq , respectively, the non-rotating e voluti on tracks 
with standard mass loss of ISchaller et al.l (l992) for stars 
with masses M < 12 Mq, and t he non-rotatin g evo lution 
tracks with high mass loss of iMevnet et ail (|l994l ) for 
stars with M > 1 2 Mq , and individual stellar spectra 
from iLeieune et al.l ll 19971) fo r stars with plane-parallel 



atmospheres, Pauldrach et al.l (|200lf > for O stars, and 



Iffillier fc Milien" (| 199Sf ) for other stars with strong winds 



The radiative and mechanical cluster feedbacks pre- 
dicted by STARBURST99 served as input to indepen- 
dent one-dimensional constant-density steady-state simu- 
lations of the emission-line spectra of metal-poor (Z = 
0.001) Hn regions, at the discrete ages 1, 2, 3, 4 and 
5 Myr. The simulations were computed with version 07.02 
of the plasm a /mole cular photoionisation code CLOUDY 
iFerland et all (|l998l ). My model cloud is a spherically sym- 
metric shell of ISM, located at radius ro, where ro is the 
outer radius of the adiabat ically expanding superbubble of 
iMac Low fc McCravl | 1988h . I let the volume of the super- 
bubble be Vb and that of the Hn region be V\. In order 
to explore why giant Hn regions do not show significant 
metal enhancements attributable to contamination with the 
ejecta from embedded massive stars, I modeled their self- 
enrichment in oxygen ejected by WR stars. This was done by 
uniformly and instantaneously enriching the volume Vb + Vi 
of valid models, with stellar ejecta predicted by STAR- 
BURST99. 

I defined a valid Hn region model as one that satisfies 
the following three conditions. First, the Stromgren radius 
of the Hn region must be larger than ro, since the Hn re- 
gion was placed at ro. Second, the time-scale to reach ther- 
mal and ionisation balance must be i cq < 1 Myr, since the 
models are 1 Myr apart. Third, the value of the [O n]/[0 ill] 
emission-line strength ratio, R23 = 3727/(4959 + 5007), 
must be in agreement with observed values in metal-poor 
giant Hn regions. 

I let the contribution of massive star ejecta to the total 
oxygen abundance of the enriched Hn region be Alog(0/H). 
The effect on Alog(0/H) of changing the dust content of the 
cloud is negligible, due to the low metallicity of the mod- 
els. The effect on Alog(0/H) of changing the value of n(H) 
by a few xlOcm -3 is also negligible. On the other hand, 
reducing the value of ro by a few xlOOpc, significantly in- 
creases the value of Alog(0/H). It would be interesting to 
include the effect of gravity in the model of the expansion 
of a superbubble, since it would reduce the value of ro and 
increase the value of Alog(0/H). My main result is that 



in metal-poor Hn regions, self-enrichment in oxygen from 
WR stars is unlikely to be the cause of observed variations 
in O/H, larger than 0.025 dex, over spatial-scales of a few 
10 2 pc. This is because for objects similar to those modeled 
in this paper, self- enrichment due to the winds of Wolf-Rayet 
(WR) stars results in a maximum oxygen abundance en- 
hancement Alog(0/H) = 0.025 dex. The latter enhancement 
is produced by cluster 1 at 4 Myr, and it is within the range 
of uncertainties i n the O/H abund ances of metal-poor giant 
Hn regions (e.g. lNava et al.ll200fj ). 

Prior to computing simulations of Hn region self- 
enrichment at ages t > 5 Myr, the following should be con- 
sidered. First, if the ionising cluster is like cluster 1, then as 
time increases, the H n region becomes fainter and harder 
to observe, but even worse, for t > 5 Myr, t cq > 1 Myr, and 
also, the value of f?23 becomes larger than observed values 
reported in the literature. In addition, at times t > 5 Myr, 
and for the two star formation laws under consideration, the 
stellar ejecta may not be confined within radius ro anymore, 
since the superbubble may have broke open into the halo of 
the host galaxy. Finally, at times t > 5 Myr, it is necessary 
to determine if the SN mechanical energy transformed into 
ionising diffuse radiation needs to be included as an addi- 
tional source of ionisation (e.g., see lCho fc Kangll2008f ). 

In the future, it would be interesting to study H II re- 
gion self-enrichment with time-dependent CLOUDY models, 
and also, to include in the study, elements such as He, C and 
N. Also, in order to complete this study, one should analyze 
the effects of changing the mass and space distribution of 
th e stars. In part i cular, one should keep in mind the work 
of lErcolano et all (|2007l ). who computed using six different 
methods the oxygen abundances that would be observed for 
spherical HIIRs of constant density, n(H) = 100 cm -3 , ra- 
dius, ro = 9 pc, metallicities Z — 0.020 and Z = 0.001, and 
ionised by stars with different spatial distributions, with a 
fixed total H-ionising photon emission rate, Qh = 3.8 x 10 50 
s _1 , such that half of the power is provided by stars with 
masses of M — 37 M q, and the othe r half by stars with 
masses of M = 56 Mm. lErcolano et al.l (|2007ft found i) that 
the typical difference in log(0/H) between the limiting cases 
where the stars occupy the same point in space, as in the 
present work, or are randomly distributed within the spher- 
ical cavity of radius ro = 9 pc, ranges from 0.1 — 0.3 dex, 
and ii) that the value of log(0/H) that would be measured, 
would be smaller for the point source scenario, for five of the 
abundance-determination methods adopted. 
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APPENDIX A: SUPERBUBBLE MODEL 

As mentioned in the introduction, massive star winds and 
SNe sweep up and compress the ambient ISM into shells. If 
the mechanical energy producing the shell is provided by one 
star, then the shell is called a bubble, and if it is provided 
by several stars, then it is called a superbubble. The size 
of a single star bubble is a fe w xlO pc, where as that of a 
sup erbubble is a f e w x 10 2 pc (|Chu et al.ll2004r i. 

IWeaver et al] (|l977l ). henceforth WMC77, have de- 
scribed the structure and evolution of an idealized bubble, 
which expands through a uniform ISM, driven by a spheri- 
cally symmetric wind with a constant mechanical luminos- 
ity. Starting from the inside, the three distinct regions of this 
bubble, rl-r3, and the boundaries of each region, bl-b3, are: 
(rl) a freely expanding wind, (bl) a wind termination shock, 
(r2) a shocked wind region, (b2) a contact discontinuity re- 
gion, (r3) a shocked ISM region, and (b3) a forward shock 
propagating into the ambient ISM (see fig. 1 in WMC77). 
In addition, the bubble's evolution can be divided into four 
stages, sl-s4: (si) free expansion of the wind for the first 
few hundred years, (s2) adiabatic expansion of the bubble 
for the next few thousand years (since the expansion is rapid, 
radiative losses do not have time to affect any part of the 
system), (s3) snowplow phase for most of the stellar lifetime 
(this phase starts when radiative losses cause the ISM swept 
up by the bubble to cool and collapse into a thin shell, and 
during this phase, the shocked stellar wind still conserves 
energy), and (s4) time when radiative losses also affect the 
dyn amics of the shocked w ind re gion . 

iMac Low fc McCravl (|l988l 'l. henceforth MM88, ex- 
tended the work of WMC77 to the multiple-star case, by 
considering the superbubble blown by the winds and SNell 
of stars in an OB association. In the model of MM88, 
the superbubble also expands through a homogeneous ISM, 
since according to the authors, the interaction of the SN 
ejecta with the pre-existing wind bubble does not affect 
the later dynamics of the superbubble. In addition, the 
superbubble's structure is similar to that of the bubble 
described in WMC77, since the stellar ejecta is assumed 
to be thermalized and co nfined within the superbubble. 
IMac Low fe McCravl (|l988h concentrate on the evolution of 
the superbubble subsequent to the collapse of the swept up 
ISM into a thin cold shell, as in stage (s3) mentioned above. 



For our parameters, this collapse occurs a few ~ 10 4 yr af- 
ter the beginning of star formation, as given by eqn. (8) in 
MM88. Throughout this paper, I considered that the super- 
bubble is in stage (s3). As in WMC77, the contact discon- 
tinuity region of the superbubble is dominated by thermal 
conduction from the hot interior to the cold shell, and the 
flux of conductive energy is balanced by the flux of mechan- 
ical energy associated with evaporation of shell gas, into 
the interior. Accordingly, the superbubble is filled with gas 
which is a combination of stellar ejecta and ISM. 



APPENDIX B: ABSORPTION WITHIN V 

For the purpose of radiative transfer, the volume of radius 
ro (Vb) was treated as a cavity empty of gas. Here I show 
that the number of H-ionising starlight photons, Q(H), is 
not greatly reduced by absorptions within Vb. 

First, I note that H is the domin ant absorber of ionising 
radiat ion within Vb- As discussed in lOsterbrock fc Ferlandl 
HqqE § 2.4), in the low-density limit, which holds within 
Vb, absorptions by He do not greatly reduce Q(H). In ad- 
dition, due to the low metallicity of the models, absorp- 
tions of ionising photons by metals do not greatly reduce 
Q{H) either. Second, I assume that each hydrogen ionisa- 
tion by a starlight photon is exactly balanced by a hydro- 
gen recombination. Therefore, showing that Q(H) is not 
greatly reduced within Vb, is equivalent to showing that 
the number of hydrogen recombinations within Vb, Ao(H), 
is much less than Q(H). Third, I use the approximations: 
n ~ n(H)+n(He) « 1.1 n(H), n e ~ ii p +u(He + ), n p ps n/1.1, 
and n(He + ) « 0.1 n, where n, n e , n p , and n(He + ) are the 
atomic, electron, proton and He + number densities, respec- 
tively. These approximations yield: n e n p « 0.92n. There- 
fore, the value of Ao(H) is given by: 

ylo(H) = 47r / a^(H, T) (0.92 n 2 ) r 2 dr, (1) 
Jo 

where is the radius where the superbubble's tempera- 
ture has dropped to 10 4 K (which is the typical equilib- 
rium temperature of an Hn region), and Q^(H, T) (cm 3 
s" 1 ) is the case A radiative recombination coefficient to form 
Hi in any energy level. For a r A , I used the parametrization 
of I Verner fc Ferlandl (|l99rj ). For n and T I used Eqns. © 
and respectively. I found log^lo(H) < 46.94 s _1 for all 
my models. On the other hand, logQ(H) > 51.95 s _1 for 
all star clusters and times considered (see Table [2}. Since 
Ao(H) <C Q(H), the gas within Vb does not significantly 
absorb the ionising radiation from the stars. 

Note that Eqn. (??) gives a lower limit to the den- 
sity inside the superbubble, since it assumes that the 
interior is dominated by t hermal conduction, whereas 
Ide Avillez fc Mac Low! (2002) strongly suggest that turbu- 
lent mixing dominates at the contact discontinuity region, 
i.e., they suggest that in reality, the density is higher and 
less homogeneous. However, the process of turbulent mixing 
is still under investigation and for now, I adopted Eqn. (??). 



